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Abstract 

We study equilibrium properties of a catalytically-activated annihilation A + A — > 
reaction taking place on a one-dimensional chain of length N (N — > oo) in which 
some segments (placed at random, with mean concentration p) possess special, catalytic 
properties. Annihilation reaction takes place, as soon as any two A particles land onto two 
vacant sites at the extremities of the catalytic segment, or when any A particle lands onto 
a vacant site on a catalytic segment while the site at the other extremity of this segment 
is already occupied by another A particle. Non-catalytic segments are inert with respect 
to reaction and here two adsorbed A particles harmlessly coexist. For both "annealed" 
and " quenched" disorder in placement of the catalytic segments, we calculate exactly the 
disorder-average pressure per site. Explicit asymptotic formulae for the particle mean 
density and the compressibility are also presented. 



PACS numbers: 82.65.+r; 64.60. Cn; 68.43.De 



1 Introduction. 



In many industrial and technological processes the design of desired chemicals requires the 
binding of chemically inactive molecules, which recombine only when some third substance - 
the catalytic substrate - is present [1, 2]. Within the two past decades much effort has been 
put in understanding of the peculiarities of such catalytically-activated reactions (CARs). On 
one hand, much progress was made in answering the question why and how specific catalytic 
substrates promote reactions between chemically inactive molecules (see, e.g. Ref.[3]). On 
the other hand, considerable theoretical knowledge was gained from an extensive study of 
a particular reaction - the CO-oxidation in the presence of metal surfaces with catalytic 
properties [4] (for a recent review se„ e.g., Ref.[5]). While the first aspect [3] sheds light 
on catalysation mechanisms and may allow the calculation of K e \ - the rate at which two 
reactants react being in the vicinity of each other and a specific catalytic substrate, the 
results of Refs.[4] show that the mere knowledge of K e \ is not sufficient. As a matter of fact, 
Refs.[4] have substantiated the emergence of an essentially different behavior as compared 
to the predictions of the classical, formal-kinetics scheme and have shown that under certain 
conditions such collective phenomena as phase transitions or the formation of bifurcation 
patterns may take place [4]. Prior to these works on catalytic systems, anomalous behavior 
was amply demonstrated in other schemes [6-8], involving reactions on contact between two 
particles at any point of the reaction volume (i.e., the "completely" catalytic sysems). It 
was realized [6-8] that the departure from the text-book, formal-kinetic predictions is due to 
many-particle effects, associated with fluctuations in the spatial distribution of the reacting 
species. This suggests that similarly to such "completely" catalytic reaction schemes, the 
behavior of the CARs may be influenced by many-particle effects. 

Apart from the many-particle effects, behavior of the CARs might be affected by the very 
structure of the catalytic substrate, which is often not well-defined geometrically, but must 
be viewed as being an assembly of mobile or localized catalytic sites or islands, whose spatial 
distribution is complex [1]. Metallic catalysts, for instance, are often disordered compact 
aggregates, the building blocks of which are imperfect crystallites with broken faces, kinks 
and steps. Usually only the steps are active in promoting the reaction and thus the effec- 
tive catalytic substrate is the geometrical pattern formed by these steps. Another example 
is furnished by porous materials with convoluted surfaces, such as, e.g., silica, alumina or 
carbons. Here the effective catalytic substrate is also only a portion of the total surface area 
because of the selective participation of different surface sites to the reaction - closed pores 
or pores with very small, bottleneck entrances are inaccessible to many reacting molecules. 
Finally, for liquid-phase catalytically-activated reactions the catalyst can consist of active 
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groups attached to polymer chains in solution. 

Such complex morphologies render the theoretical analysis difficult. As yet, only empirical 
approaches have been used to account for the impact of the geometrical complexity on the 
behavior of the CARs, based mostly on heuristic concepts of effective reaction order or on 
phenomenological generalizations of the formal- kinetic "law of mass action" (see, e.g. Refs.[l] 
and [2] for more details). In this way the parameters entering the equations describing the 
observables (say, the mean particles densities) are fixed by fits to experimental data and can 
deviate from the values prescribed by the stoichiometric relations of the reactions involved. 
The important outcome of such descriptions is that they provide an evidence of the existing 
correlations in the morphology of the chemically reactive environment. On the other hand, 
their shortcoming is that they do not explain the mechanisms underlying the anomalous 
kinetic and stationary behavior. In this regard, analytical studies of even somewhat idealized 
or simplified models, such as, for instance, the ones proposed in Refs.[4], are already highly 
desirable since such studies may provide an understanding of the effects of different factors 
on the properties of the CARs. 

In this paper we study the properties of catalytically-activated annihilation A + A — > 
reaction in a simple, one-dimensional model with random distribution of the catalyst, appro- 
priate to the just mentioned situation with the catalytically-activated reactions on polymer 
chains. More specifically, we consider here the A + A — ► reaction on a one-dimensional 
regular lattice which is brought in contact with a reservoir of A partilces. Some portion of 
the intersite intervals (thick black lines in Fig.l) on the regular lattice possesses special "cat- 
alytic" properties such that they induce an immediate reaction A + A — > 0, as soon as two 
A particles land onto two vacant sites at the extremities of the catalytic segment, or an A 
particle lands onto a vacant site while the site at the other extremity of the catalytic segment 
is already occupied by another A particle. 

We present here an exact solution of this model in two cases - a case when disorder in 
placement of the catalytic segments can be viewed as annealed, and a more complex situation 
with a quenched random distribution of the catalytic segments, and show that despite the 
apparent oversimplified nature of the model it exhibits an interesting non-trivial behavior. A 
brief account of these results has been presented in our earlier short publication [9] . We note 
finally that kinetics of A + A — ► reactions involving diffusive A particles which react upon 
encounters on randomly placed catalytic sites has been discussed already in Refs.[10, 11] and 
[12], and a rather surprising behavior has been found, especially in low-dimensional systems. 
Additionally, steady-state properties of A + A — ► reactions between immobile A particles 
with long-range reaction probabilities in systems with external particles input have been 
presented in Refs.[13] and [14] and revealed non-trivial ordering phenomena with anomalous 
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Figure 1: Onc-dimcnsional lattice of adsorption sites in contact with a reservoir. Filled circles denote 
hard-core A particles. Thick black lines denote the segments with catalytic properties, (a) denotes a 
"forbidden" particle configuration, which corresponds to immediate reaction, (b) depicts the situation 
in which two neighboring A particles may harmlessly coexist. 



input intensity dependence of the mean particle density, which agrees with early experimental 
findings [15]. For completely catalytic ID systems, kinetics of A + A — > reactions with 
immobile A particles undergoing cooperative desorption have been discussed in Refs.[16, 17] 
and [18]. Exact solutions for A + A — ► reactions in ID completely catalytic systems in 
which A particles perform conventional diffusive or subdiffusive motion have been presented 
in Refs.[19] and [20], respectively. 

This paper is structured as follows: in section 2 we define the model and introduce basic 
notations. In section 3 we focus on the case of annealed disorder and derive exact closed- 
form expressions for the pressure per site, as well as present explicit asymptotic expansions 
in powers of the activity for the mean particle density and the compressibility. In section 
4 we examine the case of quenched disorder. Here, we show that the thermodynamic limit 
result for the disorder-averaged pressure per site can be obtained very directly by noticing 
a similarity between the expressions defining the pressure in the model under study and the 
Lyapunov exponent of a product of random two-by-two matrices, obtained by Derrida and 
Hilhorst [21]. We also derive an explicit expression obeyed by the averaged logarithm of the 
partition function, which is valid for any chain's length TV, and present its large- N expansion; 
we show, in particular, that the first correction to the thermodynamic limit result for the 
disorder-averaged pressure per site is proportional to the first negative power of N. Explicit 
asymptotic expansions for the mean particle density and the compressibility are also derived. 
Finally, in section 5 we conclude with a brief summary of our results and discussion. 
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2 The model. 



Consider a one-dimensional regular lattice of unit spacing comprising N adsorption sites in 
contact with a reservoir (vapor phase) of identical, non-interacting hard-core A particles (see, 
Fig.l). The reservoir is steadily maintained at a constant pressure. 

The A particles from the vapor phase can adsorb onto vacant adsorption sites and desorb 
back to the reservoir. The occupation of the "i"-th adsorption site is described by the Boolean 
variable rij, such that 

I 1, if the "i"-th site is occupied, 
rii = < 

[ 0, otherwise. 

Suppose next that some of the segments - intervals between neighboring adsorption sites 
possess "catalytic" properties (thick black lines in Fig.l) in the sense that they induce an 
immediate reaction A + A — > 0, as soon as two A particles land onto two vacant sites at 
the extremities of the catalytic segment, or an A particle lands onto a vacant site at one 
extremety of the catalytic segment while the site at the other extremity of this segment 
is already occupied by another A particle. Two reacted A particles instantaneously leave 
the lattice (desorb back to the reservoir). Any two A particle adsorbed at extremities of a 
non-catalytic segment harmlessly coexist. 

To specify the positions of the catalytic segments, we introduce the quenched variable Q, 
so that Co = Cn = and 

| 1, if the i-t interval is catalytic, % = 1, 2, . . . , N — 1, 
[ 0, otherwise. 

Now, for a given distribution of the catalytic segments, the partition function Zjv(C) °f 
the system under study can be written as follows: 

N-l 

z N (o=j2^= ini n^-G^+O' 

{m} i=i 

where the summation ^{n 4 } extends over all possible configurations {rij}, while z denotes the 
activity, 

z = exp(/?/z), (2) 

/U being the chemical potential, which accounts for the reservoir pressure and for the particles' 
preference for adsorption. Note that Zn(() in Eq.(l) is a functional of the configuration 

C = {Q}- 

It might be instructive to remark that ^at(C) can be also thought of as a one-dimensional 
version of models describing adsorption of hard-molecules [22-29], i.e. adsorption limited by 
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the "kinetic" constraint that any two of the molecules can neither occupy the same site nor 
appear on the neighboring sites. The most celebrated examples of such models are furnished 
by the so-called "hard-squares" model [22-26], or by the " hard- hexagons" model first solved 
exactly by Baxter [28]. 

These models ehxibit phase transitions. The universal classification of phase transitions 
is known to depend on the dimensionality, the presence of further interactions and the way 
in which the lattice can be partitioned into sublattices. For bipartite lattices and interactions 
dominated by nearest-neighbor exclusion, the ordering transition is the result of competition 
between the two sublattice densities. The phase transition is thus associated with a breaking 
of the symmetry between these two sublattices. For geometrically more complex Baxter's 
hard-hexagon model, which consists of particles with the nearest-neighbor exclusion on the 
triangular lattice, the phase transition belongs to the three-state Potts model universality 
class, in accordance with the fact that the phase transition is associated with symmetry 
breaking involving three competing equivalent sublattice densities. For more discussion see 
Ref.[30]. 

In our case of the CARs on random catalytic substrates the nearest-neighbor exclusion 
constraint is introduced only locally, at some specified, randomly distributed intervals. Such 
locally frustrated models of random reaction/adsorption thus represent a natural and mean- 
ingful generalization of the well-studied exclusion models over systems with disorder. Of 
course, in this context two-dimensional situations are of most interest, but nonetheless it 
might be instructive to find examples of such models which can be solved exactly in one 
dimension. 

Our main goal here is to calculate the disorder-average pressure per site: 

p ~""S^< ln(Z " (c)) }<' (3) 

where the angle brackets with the subscript C here and henceforth denote averaging over all 
possible configurations {d}- In this case, we suppose that Q are quenched, independent, 
randomly distributed variables with distribution 

p(0=pS((-l) + (l-p)S(0- (4) 

As well, we will consider the case when the disorder in placement of the catalytic segments 
can be viewed as annealed and the mean density of the catalytic segments is equal to p; in 
this case, which requires much simplier analysis, the pressure per site is given by 

We note that such a situation can be realized in practice in case when the catalytic agents, 
modelled here as the catalytic segments, diffuse. On the other hand, an assumption of the 
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annealed disorder is often used as a meaningful "mean-field" approximation for systems with 
quenched disorder. Hence, it might be instructive to consider this case in order to check the 
behavior provided by such a mean-field approach against an exact solution in the quenched 
disorder case. 

Once Pqo are obtained, all other pertinent thermodynamic properties can be readily evalu- 
ated by differentiating with respect to the chemical potential u; in particular, the disorder- 
average mean particle density n will be given by 



while the compressibility kx obeys 



noo = J^-Poo, (6) 



U 1 dn °° (T\ 



We set out to show that for both annealed and quenched disorder cases, when Q are in- 
dependent, two-state random variables all these functions can be evaluated explicitly, in a 
closed form. We will distinguish between these two cases by assigning superscripts (ann) and 
(quen). 

To close this section, we display the results corresponding to two "regular" cases: namely, 
when p = and p = 1, which will serve us in what follows as some benchmarks. In the p = 
all sites are decoupled, and one has the trivial Langmuir adsorption results: 

P£ a ") = iln(l + z), n^) = T ^-, (8) 

and 

IT 1 4^ = \ (9) 
The "regular" case when p = 1 is a bit less trivial, but the solution can be still straightfor- 
wardly obtained. In this case, we have 

pC^> = Ilnp£±5±iY n^) = l ?i 7=> (10) 

(3 V 2 y oo l + 4z- Vl + 4z 

and 

/3 _1 4 re9) = ; ; (H) 

Note that in the p = 1 case (the completely catalytic system) the mean particle density tends 
to 1/2 as z —> oo (compared to n^ c ™' ) — > 1 behavior observed for the Langmuir case), which 
means that the adsorbent undergoes "ordering" transition and particles distribution on the 
lattice becomes periodic revealing a spontaneous symmetry breaking between two sublattices. 
In the limit z —> oo the compressibility vanishes as k^ eg ^ oc 1/ y/z compared to the Langmuir 
behavior kip an ^ ozl/z. 
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3 Annealed Disorder. 



We start our analysis of the random reaction/adsorption model considering first the situation 
in which the disorder in placement of the catalytic segments can be viewed as annealed. In 
this case, the disorder-averaged pressure per site is defined by Eq.(5) and thus has a more 
simple form than that in Eq.(3), since we have to perform averaging not of the logarithm of 
the partition function in Eq.(l) but of the partition function itself. 

Averaging of the partition function in Eq.(l) over the distribution of the catalytic segments 
can be performed very directly and yields 

N-l 

Z n = (z n {Q) =Y J Z^ ni H(l- Pni n i+I y (12) 

{m} i=i 

Since (1 — p ni nj+i) = exp(ln(l —p) m rij+i), the result in the latter equation can be thought 
of as a partition function of a one-dimensional lattice gas with nearest-neigbour repulsive 
interactions of amplitude In (1/(1 —p)). Note that here the original constraint that no two 
particle can be located simultaneously at the extremeties of the catalytic segments is replaced 
by a more tolerant condition that the particles may occupy neighboring sites anywhere, but 
the penalty of 21n(l — p) has to be paid. For any finite p < 1 this penalty can be overpassed 
by increasing the chemical potential and hence, for large z one may thus expect completely 
different behavior in the annealed and quenched disorder cases. On the other hand, for p = 1 
this penalty gets infinitely large and thus p = 1 is a special point. 

Now, to find an explicit form of Zn we proceed as follows. Let us first introduce an 
auxiliary, constrained partition function of the form 

N-2 

z' N = z N \ nN=1 = z ni n i 1 - p n * ^+0 i 1 - p ^-0 ' ( i3 ) 

{Ui} i=l 

i.e. Z' N stands for the partition function of a one-dimensional lattice gas with a nearest- 
neighbor repulsion and fixed occupation of the site i = N, tin = 1. Evidently, we have 
that 

Z N = Z N ^ + Z' N . (14) 

Next, considering two possible values of the occupation variable njv-i, i-e. rajv_i = and 
njy-i = 1, we find that Z' N can be expressed through Zjy_2 and Z' N _ 1 as 

JV-3 

Z' N = z ^ E ^ 2 n% II i 1 ~ P ni ni ^) + 
K} i=i 

N-3 

+ z 2 (l -p) ^z^ 2 "* Yl (l - V ni - p nAT- 2 ) = 

{m} i=i 

= zZ N _ 2 + z(l-p)Z' N _ 1 (15) 
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Now, recursion in Eq.(14) allows us to eliminate Z' N in Eq.(15). From Eq.(14) we have 
Z' N = Zn — Zn-i, and consequently, we find from Eq.(15) that the unconstrained partition 
function Zjy in Eq.(12) obeys the following recursion 

Z N = (l + z(l -pfjZ N -i + zpZ N _ 2 , (16) 

which is to be solved subject to evident initial conditions 

Z = l and Z 1 = l + z. (17) 

Solution of the recursion in Eq.(16) can be readily obtained by standard means, i.e. evaluating 
the generating function for Zn, Z t = X)jv=i ZNt N , and then inverting it with respect to the 
variable t, which yields 

Z N = f -N _ (1 + ^-) t -N (lg) 

zpt+{t+-t-) + zpt-(t+-t-)~ 



where 



1 + z 



(1-P)) 



t±=± 2^V( 1+ z{l ~ p) ) 2+ 4zp ~ - ' 7z P " J (19) 

Noticing next that t+ < \t-\ we find that in the annealed disorder case in the thermodynamic 
limit the disorder average pressure per site is given by 

which expression is valid for any z and p. 

Consider now the asymptotic small- z and large- z behavior of the pressure p£ nn \ the 
mean density n^™^ and the compressibility k^ nn \ Expanding P^ nn ^ in Eq.(20) into the 
Taylor series in powers of the activity z, we find that in the small- z limit P^ nn ^ follows: 

pPt nn) = z-(±+p)z+(± + 2p + p 2 )z 3 -(\ + 3p + \p 2 +p 3 )z A + 0(z 5 ). (21) 
Note that P^ nn ^ in Eq.(21) reduces to 

0P g°n) =z _ l _ z+ l _ z l_ l z 4 + 0{z ^ (22) 

and 

PPt 9) = *~\* + f - 3 " f ^ + 0(z% (23) 

for p = and p = 1, respectively. From Eq.(21) we find that in the annealed disorder case in 
the small- z limit the mean particle density is given by 



nt nn) = z 



- (l + 2p)z 2 + (l + Qp + 3p 2 )z 3 - (l + 12p + 18p 2 + 4p 3 )z 4 + C(z 5 ), (24) 
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while the compressibility obeys 

p-i k (ann) = I + - p) z - 4pz 2 + 3p^2 + 3p) z 3 - 8p(l + 4p + 2p 2 ) z 4 + 0(z 5 ). (25) 

We consider next the asymptotic behavior of P^ nn ^ in the large-z limit. We notice first 
that here p = 1 is actually a special point; that is, asymptotic large-z behavior of P^ nn ^ 
is completely different for p < 1 and p = 1 (completely catalytic systems). For p < 1 and 
z » (1 — p)~ 2 , we have that the asymptotic expansion of P^ nn ^ reads 

wL-) = m (a - p)»,) - hi - P) + * - + 



(l-p) 2 z 2(1 -p) 4 z 2 
;i + 6p + 3p 2 ) (1 + 12p + 18p 2 + 4p 3 ^ 
3(1 -p) 6 z 3 4(1 -p) 8 z 4 



while in the regular, completely catalytic case p = 1 it follows 

- pT<-> - i ln( Z ) + ^ - ^ + — ^ + O(j^). (27) 

Consequences of such a difference can be seen in a dramatically different behavior of the mean 
particle density. For p < 1 and z» (1 — p)~ 2 we find 

n (ann) = 1 1 , (j + jg) 

(l-p) 2 z (l-p) 4 z 2 
(1 + 6^ + 3^2) (i + i2 p+ l8p 2 + 4p 3 ) 5 
(l-p)6z3 + (l-p)8 z 4 + 

while in the regular case p = 1 the mean particle density is given by 

»&" 4-^ + ^-^ + 0^). ( 29 ) 

This signifies, in particular, that for p arbitrarily close but not equal to unity, the mean density 
is equal to 1 as z = oo, while for p strictly equal to unity the mean density n^ nn) = 1/2. The 
behavior of function of z for different values of p is depicted in Fig. 3. 

In a similar fashion we find that asymptotic behavior of the compressibility kx is very 
different for p < 1 and p = 1. For p < 1 and z 3> (1 — p)~ 2 , k^ nn ^ obeys 



r u(H _ 4 P | 

p (l-p) 2 z (l-p) 4 z 2 



3p(2 + 3p) 8p(l + 4p + 2p 2 ) 5 

+ (i-p)6^3 (i-^ +u ^ z h m 



while for p = 1 and z 3> 1 it follows 

1 

2zV2 1 2z ' 16z 3 / 2 256z 5 /2 1 ~Vz 7 /2 



r 1 ^ = + h + ^ - JLr. + ofi) . (3D 
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Finally, we realize that in the annealed disorder case for any fixed z the compressibility 
k^p nn ^ appears to be a non — monotonia function of p. To see this, it suffices to notice that, 
first, fcy an ' ) < k!^ e9 \ i.e. for any fixed z the value of the compressibility for p = is always 
less or equal to its value for p = 1. Second, one readily finds that in the vicinity of p = 1 the 
compressibility /c^ mra - ) obeys 

r 1 ^ = 0~ 1 4 e9) + (1+ 4 ^ )3/2 (i - p) + o((i- P ) 2 ), (32) 

i.e. for any z the value k^ &9 ^ corresponding to p = 1 is approached from above. Consequently, 
for any fixed z the compressibility k^ nn ^ is a non-monotoneous function of the mean density p 
of the catalytic segments. Behavior of the compressibility k^ nn ^ as a function of p for several 
different values of z is presented in Fig.4. 



4 Quenched Disorder. 

We turn now to the more complex situation with a quenched disorder, in which case, in order 
to define the disorder-averaged pressure, we have to perform averaging of the logarithm of the 
partition function in Eq.(l). Consequently, here we aim to determine the recursions obeyed 
byZ N (() and (ln(Zjv(C)) . 

4.1 Recursion relations for Z N (C) and (^\n(Z N (()^ . 

We proceed here along essentially the same lines as in the previous section. We introduce 
first a constrained partition function of the form 

N-2 

Z' N (Q = Z N (Q\ nN=1 =zJ2^= llni II (l-C*n*nm)(l-CiV-in,v-i), (33) 

{m} i=i 

where Z' N (() now stands for the partition function of a system with fixed set ( = {d} and 
fixed occupation of the site i = N, n^ = 1- Similarly to Eq.(14), we have that ^jv(C) obeys 

Z N (Q = Z N - 1 (Q + Z' N (Q. (34) 

Next, considering two possible values of the occupation variable njv-i, i.e. n^-i = and 
n^v_i = 1, we find that Z' N (Q) can be expressed through Z^^iC) an d Z' N _ 1 {Q as 

Z' N (Q = zZ N „ 2 (C) + 2(1 - Ov-i)^-i(C), (35) 

which parallels the result in Eq.(15). Eliminating Z' N (() in Eq.(35), we find eventually that 
the unconstrained partition function Zn(() in Eq.(l) obeys the following recursion: 

Z N (C) = (l + z(l-Cjv-i))^-i(C)+<Ar-i^- 2 (C), (36) 
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which is to be solved subject to the initial conditions in Eq.(17). 

A conventional way (see, e.g. Ref.[31,32]) to study linear random three-term recursions 
is to reduce them to random maps by introducing the Ricatti variable of the form 

Rn(C) = (37) 

In terms of this variable Eq.(36) becomes 

Rn(0 = (l + z(l - Civ-i)) + ■ with Ri{Q = Ri = l + z, (38) 

which represents a random homographic relation. Once Rn(C) 1S defined for arbitrary N, the 
partition function Zn(Q can De readily determined as the product, 

N 

Z N (0=Y[Ri(0, (39) 
i=i 

and hence, the desired disorder-average logarithm of the partition function will be obtained 

as 

N 

(hxZ N (C)) =E( ln ^(0), (40) 
f i=i c 

Before we proceed further on, some comments on the recursion in Eq.(38) are in order. 
We recall first that, by definition, each quenched random variable Q assumes only two values 
- 1 (with probability p) and (with probability 1 — p). Hence, we may formally rewrite the 
random homographic relation in Eq.(38) as 



Ri(0 



1 + z/Ri-i(C), Ci-i = 1) ( witn probability p), 
1 + z = Ri, Ci-i = 0, (with probability 1 — p). 



Note now that recursion schemes of quite a similar form have been discussed already in 
the literature in different contexts. In particular, two decades ago Derrida and Hilhorst [21] 
(see also Ref.[34] for a more general discussion) have shown that such recursions occur in the 
analysis of the Lyapunov exponent F(e) of the product of random 2x2 matrices of the form 

^A^HnQ, ;)]))«„■ («> 

where z% are independent positive random variables with a given probability distribution p(z). 
Equation (41) is related, for instance, to the disorder-average free energy of an Ising chain 
with nearest-neighbor interactions in a random magnetic field, described by the Hamiltonian 

N-i N 

1=1 i=\ 
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in which one sets J' = ln(l/^/e) and h! i = ln(l / 'y/zi) . As noticed in Ref.[21], the product in 
Eq.(41) also appears in the solution of a two-dimensional Ising model with row- wise random 
vertical interactions [35], the role of e being played by the wavenumber 6. The recurence 
scheme in Eq.(41) emerges also in such an interesting context as the problem of enumeration 
of primitive words with random errors in the locally free and braid groups [33] . Some other 
examples of physical systems in which the recursion in Eq.(41) appears can be found in [31]. 
Further on, Derrida and Hilhorst [21] have demonstrated that F{e) can be expressed as 

1 N 

F(e) = lim — V (]n%) , (43) 

1=1 

where R\ are defined through the recursion 

i^ = l + z i _i + z i _i(e 2 -l)/12;_ 1 , with R[ = l. (44) 

Moreover, they have shown that the model admits an exact solution when 

p(z) = (l-p)S(z)+ P S(z-y), (45) 

i.e. when similarly to the model under study, Z{ are independent, random two-state variables 
assuming only two values - y with probability p and with probability I— p. Supposing that 
when i increases, a stationary probability distribution P(R') of the R\ independent of i exists 
[36], Derrida and Hilhorst [21] have found the following exact result: 

F(e) = pIii(l + 6)-p(2-p)ln(l + 6^£^) + 



+ a-^E^Mi+Kf^)^ 1 )' ( 46 ) 

where 



N=l by/ 



6=1 + 



(1-y) r / e 2 y \i/2 

yj 1-1 + 4 y » 



(47) 



2e 2 y L V {l-yf 
In particular, Eq.(46) shows a striking behavior in the e — > limit; in this case, Derrida and 
Hilhorst [21] have demonstrated that for 

py > 1, and p < 1, (48) 

which implies that J p(z) ln(z) < 0, the Lyapunov exponent F(e) defined by Eq.(46) exhibits 
an anomalous, singular behavior of the form 

F(e)~e°, where a = - ln(p)/ln(y). (49) 

We turn now back to our recursion scheme in Eq.(38) and notice that setting 

R i (0 = (l + z)R' i , (50) 
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and choosing 

y = - T ?L- = -n«r ) , and e 2 = = n£ n \ (51) 

makes the recursion schemes in Eqs.(38) and (44) identic! Consequently, the disorder-average 
pressure per site in our random catalytic reaction/adsorption model can be expressed as 

p^en) = Iln(l + z) + ±F(e), (52) 

where F(e) is the Lyapunov exponent of the product of random 2x2 matrices in Eq.(41), in 
which e and Z{ are defined by Eqs.(45) and (51). 

Note next that the first term on the right-hand-side of Eq.(52) is a trivial Langmuir re- 
sult for the p = case (adsorption without reaction) which would entail = z/(l + z). 
Hence, all non-trivial, disorder-induced behavior is embodied in the Lyapunov exponent F(e). 
We hasten to remark, however, that despite some coincidence of results, the random reac- 
tion/adsorption model under study has completely different underlying physics, as compared 
to the model studied by Derrida and Hilhorst [21]. Thus, one would not expect any singular 
overall behavior of pressure in the e — > limit (which corresponds here to the limit of vanish- 
ingly small activities z (or \i — > — oo), and thus pertains to n <C 1). In consequence, here y is 
also dependent on z and y — ► in the same manner as e. Moreover, in our case y < 0, which 
invalidates the condition in Eq.(48). 

4.2 Disorder-averaged pressure. 

The disorder-averaged pressure per site can be thus readily obtained from Eqs.(46) and (47) 
by defining the parameters y and e as prescribed in Eq.(51). This yields the following explicit 
representation 



(3Pt en) = m(^)-(l-p)ln(l-u, 2 ) + 



a\2 00 

- ^^L^^ln(l-(-l)V+ 2 ), (53) 



P N=l 



where 



and 



1 + y/l + 4z 



(54) 



VT+^-i 2 l 

lo = , = z/(pi = 1 — (55) 

v / TT4l+i <f> z 

Note that <f> z obeys (pzi^z — 1) = z and thus for z = 1 the <t>\ is just the "golden mean", 
<pi = (\/5 + l)/2. Below we will show why and how this mathematical constant appears here. 

On the other hand, the derivation of the result in Eq.(53) can be performed in a very 
straightforward manner without resorting to the assumption on existence of a stationary 
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probability distribution P(R'); the intermediate steps of such a derivation contain useful 
formulae, which might be helpful for the understanding of the asymptotic behavior of Eq.(53). 
Since it allows us to answer also the question of how the thermodynamic limit is achieved, 
we find it expedient to present such a derivation here. 

We start with calculation of an explicit form of < In Ri(() >£. To do it, it suffices to 
notice the following two points: First, we notice that 

^(0=^(6-1, 0-2,0-3, ...,&), (56) 

and 

^i-l(C) = Ri(Ci-2, 0-3, 0-4, • • • , Cl), (57) 

i.e. Ri_k(C) depends only on 0-fc-n with n = 1,2, . . . ,i — k — 1 and is independent of 0-ifc- 
Second, with probability 1—p the Ricatti variable is set equal to 1 + 2, i.e. to its initial value 
Ri, which is a non-random function. These two observations allow us to work out an explicit 
formula for < ln(i?j(£)) >^ which is valid for any i. 

Taking the logarithm of both sides of Eq.(38) and averaging it with respect to the distri- 
bution of random variables Q, we have: 

( ln(^(C))) c = ( In (l + z(l - C,_i) + z|g(C)) ) c - (58) 

We note next that since i?j_i(£) is independent of 0-1, we can straightforwardly average the 
right-hand-side of Eq.(58) with respect to 0-1, i- e - 

ln(i? i (C))) c = (ln(l + z(l-Ci-i)+^(C))) f = 
= (1 - p) ln(l + z) + p( In (l + ^ ) 
= (l-p)ln(l + z) + 



+ Kin 1 + ) (59) 

1+2(1 -°- 2)+2 ^(o 

Now, since that i?j_2(C) is independent of 0-2, we can again average over states of this 
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variable, which yields 



(1 - p) ln(l + z) + p( In (l + 



1 + z(l - Q- 2 ) + z 



Ci-2 J I c 



= (l-p)ln(l + 2;)+p(l-p)lii(l + Y^— ) + 

-^-) 



i*i-2(C) 



+ / In 1 + 



^- 2 (C) 

= (l-p)ln(l + z)+p(l-p)ln(l + r ^-) + 
+ p 2 (|ln(l + - 



+ 



(60) 



1 + z(l - Ci- 3 ) + 



Ci-3 



^i-3(C) 

Noticing again that 3(C) is independent of £i-3 and etc., we arrive eventually at the 



(61) 



following explicit representation for ( ln(i?j(£))y : 



i-1 



ln(i2i(C))) = (1 -P) J> n_1 .Fn +P % - l H 

^ 71=1 



where the sum on the right-hand-side of Eq.(61) is defined for i > 2 and equals zero otherwise, 
while T n denote natural logarithms of the Stieltjes-type continued fractions of the form 

/ \ 

^i=ln(l + z), T 2 = Infl + — ^-V f 3 = ln 1 + - 

\ 1 + z / I 1 



+ 



l + zj 



( 



Ti = In 



1 + 



1 + 



1 + 



V 



1 + 



(62) 



l + z 



J 



To analyze the leading large- iV behavior of the disorder-average pressure per site we resort 
to the standard generating function technique [37], often used, in particular, in the analysis of 
peculiar properties of different random walks [38]. Let us define first an auxiliary generating 
function 

00 

n t = J2t n (HRn(0)) ■ (63) 

n=l C 

Then, multiplying both sides of Eq.(61) by t n and performing summation, we readily find 
that 

(64) 



P(l-t)^ 
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Consequently, the generating function of the averaged logarithm of the partition function 
Z N (Q obeys: 

oo oo N 

Z t = ^^<lnZ JV (C)>C=E^E( ln (^(0)) = 
N=l N=l n=l ^ 

and hence, the generating function of an average pressure per site, defined as 

1 00 t N 

Vt= 3^2 AT < ln ^(0>C> (66) 

P N=l 



attains the form 



where 



^ 00 

V * = P Njr " i 1 " ~ P^+i) > (67) 



Is = I dr JT — f . (68) 

Now, in the large-iV limit, the asymptotic behavior of the disorder-average pressure P^v 
per site in a finite chain of length N can be obtained very directly from the expansion of Vt 
in the vicinity of the closest to the origin singular point [37], i.e. t = 1. Since, in the limit 
t — > l - , In obeys 

J Ar = -l_ + (iV-l)ln(l-i) + 0(l), (69) 
we have that in this limit Vt is given by 

v t = Yzr t p ™ uen) + ln (! - *) {p^ p t en) ) +0(1), (70) 

where 

p(quen) 



n=l 



Consequently, we find that in the large- N limit p^ uen ^ follows 

p (quen) _ p ( qU en) _ J_ ( ®_ p(quen)\ , /nf_L\ (70^ 

in which equation P^ uen ^ in Eq.(71) is the desired thermodynamic limit result for the disorder- 
average pressure per site in the quenched disorder case. Note that in virtue of the expansion 
in Eq.(72), the corrections to the thermodynamic limit are proportional to the first inverse 
power of the chain length N. Note also that since 



lim T n = H<f> z ) = In ( 1 + V ! + 4 ^ (73) 
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i.e. J~ n is the 7i-th approximant of ln(0 z ), i"oo 

can be thought of as the generating function 
of such approximants. One expects then that for z < 1 the sequence of approximants converges 
quickly to ln(0 z ); expanding the n-th approximant T n into the Taylor series in powers of z, 
one has that the first n terms of such an expansion coincide with the first n terms of the 
expansion of ln(4> z ). Consequently, T n and T n -\ differ only by terms of order z n , which 
signifies that convergence is good. On the other hand, for z > 1 convergence becomes poor 
and one has to seek for a more suitable representation. As a matter of fact, already for z = 1 
one has that in the limit n — > oo the approximant T n tends to ln(0i), i.e. the logarithm of the 
"golden mean", which is known as the irrational number worst approximated by rationals. 
Moreover, for z — ► oo the convergence is irregular in the sense that only the approximants 
with odd numbers show the same large- z behavior as ln(4> z ); the approximants with even n 
all tend as z — > oo to finite values ln(n/2 + 1) (see, Fig. 2). We turn now back to the result 



3.5 -i 




5~ 10 15 20 25 30 



Z 

Figure 2: Plot of the approximants T n , n= 1, 2, 3, 4, 5 and 6, and ln(0 z ) versus activity z. 

in Eq.(71) aiming to find a convenient representation more amenable for further analysis. To 
do this, let us note that T n in Eq.(62) can be expressed as the logarithm of the convergents 
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of the Stiltjes-type continued fractions: 

where K n (z) are polynomials of the activity z defined through the three-term recursion 1 

K n {z) = K n _ 1 {z) + zK n _ 2 {z), K (z) = l, K 1 (z) = l + z. (75) 

These polynomials can be, of course, obtained very directly by introducing their generating 
function, but we can avoid doing it by merely noticing that they are simply related, in view of 
the form of the recursion in Eq.(75), to the so-called golden or Fibonacci polynomials F n+ 2{x) 
[39], which are defined by the three-term recursion of the form 

F n+1 (x) = xF n (x) + F n _ 1 (x), iq(x) = l, F 2 (x) = x. (76) 

On comparing the recursions in Eqs.(75) and (76), one infers that 

K n (z) = z(« +1 )/ 2 F n+2 (l/v^). (77) 

Hence, the approximant T n can be expressed as 

^ = 2 ln(z) + ln (?— (T77i)> (78) 

Note that even at this stage one may understand where from such functions as 4> z appear in 
the expression for the disorder-average pressure in Eq.(53) (first term on the rhs). The point 
is that, similarly to the Fibonacci numbers F n = F n (l), which obey lim n ^ 00 F n jF n -\ = </>i = 
(V5+ l)/2, the ratio of two consequitive golden polynomials F n (l/^/z) and F ri _i(l/y / z) also 
converges as n — > oo to a finite limit given by the function 4> z /^/z. One expects hence that the 
rest of terms on the rhs of Eq.(53) stems from the finite-n effects and describes the relaxation 
of the logarithm of F n (l/^/i)/F„_i(l/^/i) to ln(0 z ) . 

To determine the relaxation terms, one uses the standard definition for the Fibonacci 
polynomials: 

^)^[(^)"-<-)"(^)"] 

In virtue of this formula, one finds that the ratio of two consequitive golden polynomials obeys 



J It is straightforward to check that the polynomial K„(z) is just the partition function in Eq.(l) for a chain 
of length n in the completely catalytic p = 1 system, i.e. K n (z) = Z n (( = 1). 
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where uj has been defined in Eq.(55). Consequently, we find that the n-th approximant T n is 
given by 

1 - (-l) n io n+2 ^ 



T n = ln(cf) z ) + In 



(81) 



1 + (-l) n to n+1 

where, as we have already remarked, the first term on the rhs of Eq.(81) corresponds to 
the limiting form of the approximants, while the second term determines the relaxation to 
this limiting form; more specifically, to the leading order this relaxation is described by an 
exponential function exp(— nln(l/u;)). Consequently, one expects a fast convergence in case 
when z is small [oj is small) and poor convergence when z — > oo [oj — > 1). Substituting 
Eq.(81) into Eq.(71) we recover, upon some straightforward algebra, the result in Eq.(53). 

4.3 Asymptotic behavior of the disorder-average pressure, mean density 
and the compressibility. 

Consider first the small-z behavior of the disorder-average pressure per site, defined by 
Eq.(53). As we have already remarked, expanding the n-th approximant J- n into the Taylor 
series in powers of z, one has that the first n terms of such an expansion coincide with the 
first n terms of the expansion 

MM = I- ( 1 + v/rT ^) =~t ( ~ 1>,T( " + ! /2> («) 
yYJ V 2 J 2J^^ Tn + 1) n V 1 

v n=l 

which implies that T n and T n -\ differ only by terms of order z n and allows to obtain very 
directly a convergent small- z expansion of the pressure P^ uen ^ . We find then 



ppgien) = z _(l +p J z 2 + ^ + 2p + p 2y 

- (\ + Ip + ±P 2 +P 3 )z 4 + 0(z 5 ). (83) 
Consequently, in the small-z limit the mean density obeys 

n (g«en) = z _ Q + 2p ) z 2 + (l + Qp + 3 p A / - (\ + Up + 16p 2 + 4p 2 ) Z 4 + 0( Z 5 ), (84) 



while the compressibility k!^ uen ^ follows 

p-i k (wen) = ^ +p ^-p)z-4p(2-p)z 2 + 3p(8-p-2p 2 ^z 3 + 0(z 4 ). (85) 

Note now that the expressions in Eqs.(83) to (85) differ from their counterparts obtained in the 
annealed disorder case, Eqs.(21), (24) and (25), only starting from the terms proportional to 
the fourth power of the activity z. On the other hand, the coefficients in the small-z expansion 
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nonetheless coincide with the coefficients in the expansions of and P^ e9 ^ when we set 

in Eq.(83) p = or p = 1. 

Now, we turn to the analysis of the large- z behavior which is a bit more complex than 
the z <C 1 case and requires understanding of the asymptotic behavior of the sum 

oo 

S= 5> iV ln(l-(-l) J V v+2 ) (86) 

N=l 

entering Eq.(71). We note first that in this sum the behavior of the terms with odd and even 
N is quite different and we have to consider it separately. 
Let 

^ oo 

S odd = - Y,P 2N ^{l+" 2N+1 ) (87) 

P N=l 

denote the contribution of the terms with odd N. Note that when z — > oo (i.e. uj — > 1) the 
sum S dd tends to pln(2)/(l — p 2 ). The corrections to this limiting behavior can be defined 
as follows. Expanding the logarithm In (l+UJ 2N +^ into the Taylor series in powers of uj and 
then, using the definition to = 1 — l/(f> z and the binomial expansion, we construct a series in 
the inverse powers of <j> z . This yields 

S odd - ln(2)- - ^_^- + 

1 p(3 + 6p 2 - p 4 ) 1 1 p(15 + 10p 2 -p 4 ) 1 / J_ 
+ 8 (l-p 2 ) 3 0| + 24 (1-p 2 ) 3 0| + V0| 

Note that this expansion is only meaningful when (f> z S> (1 — p) -1 , (z S> (1 — p)~ 2 ), which 
signifies that p = 1 is also a special point for the quenched disorder case. 



Further on, plugging into the latter expansion the definition of <j> z , 4> z = (1 + y/1 + 4z)/2, 
we obtain the following expansion in the inverse powers of the activity z: 

^ " I" ^2 ln ^) -2(l-p2)2 ^2 + 

p (9-2p 2 +p 4 ) l p (3-4p 2 +p 4 ) 1 / 1 \ 

+ 8 (1-p 2 ) 3 * 48 (1-p 2 ) 3 z 3 / 2 + \z 2 )' [ j 

Consider next the sum 



oo 



^ n =^^ln(l-^+ 2 ), (90) 
JV=1 

which represents the contribution of terms with even N. Note that in contrast to the behavior 
°f S dd, the sum in Eq.(90) diverges when z — > oo (uj — > 1). Since 1 — w 2Ar+2 ~ 1 — cj for 
a; — > 1, we have that in this limit 

p 2 

Seven ~ " ~ m(l - (91) 

1 — J) z 
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Figure 3: The mean density of adsorbed particles versus the chemical potential /3^t for the annealed 
(curves tending to unity) and quenched disorder case for different values of the mean density p of the 
catalytic segments. 



To obtain several correction terms we make use of one of Gessel's expansions [40]: 

v ; fc=i 

where gk(2N + 2) are the Dedekind-type sums of the form 

g k (2N + 2)= Yl 1 . fe . ( 93 ) 

C 2Ar+ 2=1)f#1 ^ - lj 

where the summation extends over all ( being the (2iV + 2)-th roots of unity (with ( = 1 
excluded). As shown in Ref.[40], the weights gk(2N + 2) are polynomials in N of degree at 
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Figure 4: The compressibility /? _1 A;t versus the mean density p of the catalytic segments for several 
values of the activity z, z — 0.2, 2 and z — 20. Upper, non-monotoneous curves show the behavior of 
(3~ x kT in the annealed disorder case, while the lower curves correspond to the solution in the quenched 
disorder case. 

most k with rational coefficients; first few values of gu{2N + 2) are: 



Now, setting x = l/(f) z in the expansion in Eq.(92), plugging it to Eq.(90) and performing 
summations over N, we find that S even can be written down as 



5l (2iV + 2) 
g 2 (2N + 2) 
g 3 (2N + 2) 
g 4 (2N + 2) 



-(2iV + l)/2, 

-(2JV + l)(2JV-3)/12, 

(2JV + l)(2JV-l)/8, 

(27V + l)(8iV 3 + 28iV 2 - 186iV + 45)/720. 



(94) 





CO 




s, 



'even — 



2 ln (^) + 



-ln(2) +s p -Y,G k (p) 



k4> k z 



(95) 



k=i 
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where s p is an infinite series of the form 2 



s P =Y,P 2N HN + l), (97) 

JV=1 

while Gk{p) are the generating functions of the polynomials gk(2N + 2): 

oo 

Gk(p)= Y,9k(2N + 2)p 2N (98) 

JV=1 

Inserting next the definition of (f> z , we find the following explicit asymptotic expansion 



1 p 2 p 2 

Seven = ~ 7, J3^2 ln ( Z ) + ^ + Sp + 

p 2 (2 - p 2 ) 1 p 2 (21 - 18j3 2 + 5p 4 ) 1 p 2 (2 - p 2 ) 1 

"I" ?TT7^ o\ q I" TTTT^ oTo o /o "I - ^ 



(1). (99) 



(l_ p 2)2 z l/2 24(1 -p 2 ) 3 Z 24(1 - p 2 ) 2 ^3/2 

Finally, combining the expansions in Eqs.(71), (89) and (99), we find the desired large- z 
expansion for the disorder-averaged pressure P^ uen ^: 

Note that P^ en ^ in Eq.(lOO) shows a completely different behavior compared to its coun- 
terpart in the annealed disorder case already in the leading term in the large- z expansion. 
Note also that here p = 1 appears to be a special point and thus the expansion in Eq.(lOO) 
becomes meaningless for p = 1. As a matter of fact, for p arbitrarily close but less than unity 
one has intervals which are devoid of the catalytic segments. Contribution of such intervals 
to the overall disorder-average pressure is of a Langmuir-type and vanishes only when p is 
strictly equal to unity, which implies that also here p = 1 is a special point. 
We find next that for z S> (1 — p)~ 2 the mean particle density obeys 

n (quen) = __ _ _ _ _±_P_ P_ - + O (-) , (101) 

i.e. contrary to the behavior of the mean particle density in the annealed disorder case, 
Eq.(28), n^ en ^ tends towards a constant value 1/(1 +p), which depends on p and coincides 
2 Note that s p shows a non-analytic behavior when p — ► 1. This function can be represented as 

; l n (l_ p 2 )_ P-^iZ^^n,!), (96) 



- p 2 1 — p 



where $(p 2 , n, 1) are the Lerch transcedents, <J>(p 2 , n, 1) = X^o(l + n p 2L ■ It is straightforward to find then 
that s p = — 1 _ 1 ln(l — p 2 ) — yrr + C(ln(p))i where 7 is the Euler constant. 
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with the corresponding values n^ lan ^ = 1 and n^^ = 1/2 for p = and p = 1. Behavior 
of the mean density versus the chemical potential a for the annealed and quenched disorder 
cases is presented in Fig. 3. 

Finally, from Eq.(lOl) we find that the compressibility fcy" en ^ admits the following form: 

^ 2 



6 (1 +p)(l -p 2 ) z 36 (1 +p) 2 (l -p 2 ) 2 z 2 \z 3 
which expansion also holds in the asymptotic limit z> (1 — p)~ 2 



5 Conclusions. 

To conclude, in this paper we have presented an exact solution of a random reaction/adsorption 
model, appropriate to the situations with the catalytically-activated reactions on polymer 
chains containing randomly placed catalytist. More specifically, we have considered here the 
A + A — > reaction on a one-dimensional regular lattice which is brought in contact with a 
reservoir of A partilces. Some portion of the intersite intervals on the regular lattice was sup- 
posed to possess special "catalytic" properties such that they induce an immediate reaction 
A + A — > 0, as soon as two A particles land onto two vacant sites at the extremities of the cat- 
alytic segment, or an A particle lands onto a vacant site while the site at the other extremity 
of the catalytic segment is already occupied by another A particle. For two different cases; 
namely, when disorder in placement of the catalytic segments can be viewed as annealed, and 
a more complex situation with a quenched random distribution of the catalytic segments, we 
have determined exactly the disorder-averaged pressure per site. For the annealed disorder 
case such a pressure has been found in a closed form and explicit asymptotic expansions 
in powers of the activity for the mean particle density and the compressibility have been 
obtained. In the case of quenched disorder we have shown that the thermodynamic limit 
result for the disorder-averaged pressure per site can be obtained very directly by noticing 
a similarity between the expressions defining the pressure in the model under study and the 
Lyapunov exponent of a product of random two-by-two matrices, obtained by Derrida and 
Hilhorst [21]. We have also derived an explicit expression obeyed by the averaged logarithm of 
the partition function, which is valid for any chain's length N. From this expression we have 
constructed the large-iV expansion and have shown, in particular, that the first correction 
to the thermodynamic limit result for the disorder-averaged pressure per site is proportional 
to the first negative power of N. The leading term in this expansion coincides with the one 
found from the analysis by Derrida and Hilhorst. Explicit asymptotic expansions for the mean 
particle density and the compressibility were also derived. We have demonstrated that for low 
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activities in the annealed and quenched disorder cases the coefficients in the corresponding 
expansions of the pertinent parameters in the Taylor series in powers of z coincide up to the 
order z 3 and slightly deviate from each other starting from the fourth order. On the other 
hand, expansions in inverse powers of z (large-,? behavior) are different already in the leading 
order. Most spectacular difference between the annealed and quenched disorder case have 
been observed in the behavior of the compressibility: in the annealed disorder case it appears 
to be a non-monotoneous function of the mean density p of the catalytic segments, while in 
the quenched disorder case it is a monotoneously increasing function of p. 
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